Infrared spectroscopy of small-molecule endofullerenes 
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Hydrogen is one of tire few molecules which has been incarcerated in the molecular cage of 
Ceo and forms endohedral supramolecular complex H2@C6o- In this confinement hydrogen 
acquires new properties. Its translational motion becomes quantized and is correlated with its 
rotations. We applied infrared spectroscopy to study the dynamics of hydrogen isotopologs 
H2, D2 and HD incarcerated in Cgo- The translational and rotational modes appear as side 
bands to the hydrogen vibrational mode in the mid infrared part of the absorption spectrum. 
Because of the large mass difference of hydrogen and Cgo and the high symmetry of Cgo the 
problem is identical to a problem of a vibrating rotor moving in a three-dimensional spherical 
potential. The translational motion within the Cgo cavity breaks the inversion symmetry 
and induces optical activity of II2. We derive potential, rotational, vibrational and dipole 
moment parameters from the analysis of the infrared absorption spectra. Our results were 
used to derive the parameters of a pairwise additive five-dimensional potential energy surface 
for Il2@C6o. The same parameters were used to predict H2 energies inside C7o[Xu et al., 
J. Chem. Phys., 130, 224306 (2009)]. We compare the predicted energies and the low 
temperature infrared absorption spectra of H2@C7o. 



INTRODUCTION 

A small cavity inside the fullerene cage is a potential trapping site of atoms and has attracted 
attention of scientists from the moment of discovery of Ceo [l] • The demonstration of formation of 
La@C6o after laser bombardment of La-impregnated graphite was immediate [3]. Since then the 
field of studies of endohedral fullerenes has expanded. Endohedral fullerenes with noble gas (He 
and Ne in [3], Ar, Kr, Xe j3]), nitrogen [s^] or phosphorus atoms and with metal clusters 0] are 
made under extreme conditions using arc discharge, ion bombardment or high pressure and high 
temperature treatment. 

Extreme methods are not suitable for encapsulation of small molecules. A different approach, 
"chemical surgery" , was applied by Rubin when he made first open-cage fullerene ^] with a orifice 
large enough to load it with ^He or H2 using less extreme temperature and pressure [9j. Soon Y. Mu- 
rata, M. Murata and K. Komatsu synthesized another opened-cage derivative of Cgo and achieved 
100% yield in filling with H2 [l3]- There a generation of closed-cage H2@C6o was observed in the 
process of matrix-assisted laser desorption/ionization time-of- flight mass spectrometry analysis of 
this opened-cage complex. Chemical methods were developed to close opened-cage fullerenes and 
H2@C6o was produced in milligram quantities [O, To accommodate two hydrogen molecules 
a cavity larger than Cgo is needed. Two H2 were trapped in opened-cage C70 with the yield 3:97 
in the favour of species with one H2 per cage [13]- The restoration of closed cage retains approx- 



imately the same ratio of (H2)2@C7o to H^Q Crn |l4[ |. Molecules other than hydrogen trapped in 
opened-cage fullerenes are carbon monoxide 15], water 0, [13]) ammonia 18] and methane 19]. 
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Recently Kurotobi and Murata succeeded in closing one of them and making the first closed-cage 
endohedral complex with a trapped polar molecule, H20@C6o [20^- The rotational modes of endo- 
hedral water were observed by inelastic neutron scattering (INS), far infrared (far-IR) and nuclear 
magnetic resonance (NMR) at cryogenic temperatures [2ll |. 

Up to date H2@C6o has been the most studied small molecule endofuUerene. The local envi- 
ronment of a hydrogen molecule in the fullerene cage has a negligible inhomogeneous distribution 
of interaction parameters. All the trapping sites are similar in H2@C60i except the isotopic distri- 
bution of carbon atoms and crystal field effects in solid H2@C6o- The H2@Cgr) is a stable complex 
and can survive a short period of heating up to 500 °C under vacuum [12]. These properties 
make H2@C6o appealing for spectroscopic and theoretical investigations of interactions between 
the molecular hydrogen and carbon nanosurfaces. 

Three spectroscopic techniques NMR, INS and infrared (IR) , have been used to study endohedral 
hydrogen. NMR studies cover spin lattice relaxation rates of H2@C6o in organic solvents 221. |23||. 
and in the presence of paramagnetic relaxants 24 , [2H] . NMR was used to follow the ortho-para 
conversion in H2@C6o in the presence of molecular oxygen at 77 K [2^ or upon photoexcitation of a 
C70 triplet state |27|. NMR study of micro-crystalline H2@C6o samples at cryogenic temperatures 
shows splitting of the J = 1 rotational state 28, |2^, a sign of the symmetry reduction from the 
icosahedral symmetry in the solid phase. Similarly, splitting of the ground ortho state was deduced 
from the heat capacity measurements [3o|. 

The overview of the low temperature NMR, INS and IR work on H2@C6o is given by Mamone 
et al. 311] • The first IR study of H2@C6o was limited to 6K 3j]. The translational and rotational 
transitions appeared as sidebands to the hydrogen molecule bond-stretching vibrational transition, 
= — 7- 1, in the mid-IR spectral range. The direct translational and rotational transitions were 
not observed in the far-IR below 200 cm~^ 33|. The extension of IR studies to higher temperature 
made possible to probe the hydrogen-Ceo interaction potential in the ground v = and first excited 
V = 1 vibrational states and a whole range of hydrogen isotopologs H2, D2 and HD was studied 
H, 34]. The isotope effects and translation-rotation coupling were also studied by INS in H2@C6o 



and HD@C6o 35|]. The translational and rotational energies of H2@C6o and HDQCgo in the v = 
state determined by IR spectroscopy are consistent with the low temperature INS results |35l] . 
There is no Raman data on H2@C605 except a report on H2 inside an opened-cage fullerene |36l |. 

In this paper we will review the IR studies of hydrogen isotopologs in Cgo and present the 
analysis of IR low temperature spectra of H2@C7o. The far-IR properties of H20@C6o will not be 
reviewed here [21|. 



THEORY 

Quantum statistics plays an important role in the dihydrogen wavefunction symmetry and has 
a pronounced effect on the rotation of the hydrogen molecule [37f]. The symmetry relative to the 
interchange of two protons dictates that there are two forms of molecular hydrogen, called para- 
and ortho-]i2. The two proton spins {Ip = 1/2) are in the anti-symmetric / = total nuclear 
spin state in para-ii2 and in the symmetric 1 = 1 state in ortho-li2- Even rotational quantum 
numbers J are allowed for para-Ii2 and odd J for ortho-Ii2- The nucleus of D is a boson, nuclear 
spin Id = 1. Thus the rotational state with an even quantum number J has D2 in the state where 
the total nuclear spin of D2 is either zero or two, / = 0, 2. This is called ortho-T>2, while para-D2 
has the total nuclear spin 1 = 1 and odd J values. 

A homonuclear diatomic molecule with the total nuclear spin 1=1 that is rotating in its ground 
state. This J = 1 rotational state is 118 cm~^ for H2, and 58 cm^^ for D2, above the rotational 
ground state J = of even-/ species. A thermal transition J = 1 — )■ must be accompanied by 
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the change of the total nuclear spin of the molecule, a process of very low probability. The time 
constant of thermal relaxation between the ortho and the para manifolds is very long and the room 
T ortho-para ratio is maintained even at cryogenic temperatures. The equilibrium distribution of 
H2 nuclear spin isomers is rio/rip = 3 and of D2 is rio/rip = 2 at room temperature. To change 
the total nuclear spin of a molecule the two nuclei must experience different magnetic field. The 
ortho-para conversion can be activated by using a paramagnetic center as a source of the magnetic 
field gradient. The equilibrium rio/rip = 1 is reached at 77 K by dispersing H2@C6o on a zeolite 
surface and exposing it to molecular oxygen which acts like a spin catalyst [38]. There are no ortho 
and para species for HD. All rotational levels of HD are in thermal equilibrium and there is one 
rotational ground state, J = 0. 



Quantum chemistry calculations are challenging for a hydrogen molecule in a weak van der Waals 
interaction with a large fullerene molecule. The availability of experimental data on endohedral 
H2 has stimulated theoretical work in this direction. Theoretical investigations currently cover 
the calculations of rotation-translation energies of hydrogen isotopologs in Cgo fsol, Eol] and H2 in 



C70 [4l|, and the stability of Cqq or C70 with one or more incarcerated H2 142, |43|]. Empirical 



parameters of the Morse potential between H-H and contact Dirac interaction between H-C were 
adjusted 44] and density-fitting local M0ller-Plesset theory tested [iH] using the experimental H2 
vibrational frequency inside Cqq. Classical molecular dynamics and density- functional theory have 
been combined to reproduce accurately the NMR chemical shift of in H2@C6o [461] . 



Both the hydrogen molecule and fullerene have closed shell electronic structures and therefore 
the interaction between them is the van der Waals interaction. The simplicity makes H2@C6o ideal 
for the studies of non-covalent bondings between H2 and carbon nano-surfaces, the knowledge 
needed for the design of carbon-based hydrogen storage materials. The high icosahedral symmetry 
of the Ceo cavity is close to spherical and therefore H2@C6o represents a textbook example of a 
body moving in a spherical potential well 47, |4^. In addition, H2 rotates around its center of 
mass. H2 is not spherical and therefore its interaction with the walls of the cavity depends on 
its orientation what leads to the coupling between translational and rotational motion [49i] . If the 
translational and the rotational motions are coupled then in the spherical potential the conserved 
angular momentum is the sum of translational and rotational angular momenta [sol. [5o|. H2@C6o 
is a rare example where the quantum dynamics of a diatomic rotor in a confined environment 
can be studied. Another, but with limited degrees of freedom, interestin g ex ample related to the 



fullerenes is the quantum rotor C2 in a metallofullerene C2@Sc2C84 [5l|, [S^]. The two scandium 
atoms limit the translational motion and fix the rotational axis of C2 relative to the fullerene cage. 
At low temperature the rotation of C2 is hindered because it has a small rotational constant and 
is therefore more susceptible to the corrugations of the carbon surface. H2 provides examples of 
two-dimensional rotors, like H2 on a Cu surface [s^] or II2 in intercalated graphite 54]. 



High pressure loading of solid Cfin creates interstitial H2 . Exohedral II2 has been studied by IR 



0, a, INS 0, a, NMR 0, a and Raman |H spectroscopies. Hydrogen is trapped m an 
interstitial site of octahedral symmetry and theory predicts translation-rotation coupling 
However, broadening of experimental lines has prohibited accurate determination of the H2-C60 
interaction potential. 



The observed IR spectra of hydrogen encapsulated in Ceo consist of several absorption lines. 
We construct a model Hamiltonian and a dipole moment operator with few adjustable parameters 
to describe accurately the position and intensity of such multi-line spectrum. 
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FIG. 1. Coordinates of a diatomic rotor inside the Cgo cage are the relative position s of two nuclei in the 
diatomic molecule and the diatomic center of mass displacement R from the Ceo cage center. In case of 
heteronuclear diatomic molecule (HD) the center of mass is not in the middle of the chemical bond between 
the two nuclei. Two polar spherical coordinates, Q and 9, are shown. (Online version in Colour) 



Diatomic molecule in a spherical potential well 

To describe the motion of a hydrogen molecule inside Cgo we use the following model. The Cgo 
is considered to be rigid, its center of mass is not moving and it does not rotate. We treat H2@C6o 
as an isolated molecular complex and approximate the true icosahedral symmetry of an isolated 
Ceo with spherical symmetry. It means that in this approximation H2 moves in a rigid spherically 
symmetric bounding potential provided by the cavity of Cgo- Beside the translational movement 
inside Cgo the hydrogen molecule has its internal degrees of freedom, vibration and rotation of two 
nuclei relative to its center of mass. There are no coupling terms between ortho- and para-states 
in our model Hamiltonian. 

The theoretical work of Olthof et al. (63| is a comprehensive description of the dynamics of a 
loosely bound molecule inside Cgo- Olthof et al. model the intermolecular potential as a sum of 
atom-atom potentials and expand it in spherical harmonics. They determined the radial part of the 
wavefunction with discrete variable representation method. The radial part of the wavefunction in 
our approach is given by algebraic functions, solution of the three-dimensional spherical oscillator 



47l . |48| . |64| |. The advantage is that matrix elements are calculated in algebraic form avoiding 
time-consuming numerical integration. 

The position and orientation of the H2 molecule is given by spherical coordinates R = {R, 0}, 
Q = {0, and s = {s,r2s}, = {(^,4'} where R is the vector from the center of the Ceo cage 
to the center of mass of H2 and s is the internuclear H-H vector, as shown in Fig.[T] The center 
of mass translational motion of H2 is given by eigenfunctions of the isotropic three-dimensional 



harmonic oscillator j48|, l64 1 



^-^LmM^) = rNliR) YLM^m, (1) 

where V'tvl is the radial wave function and Ylml is the spherical harmonic. The size of the 
H2 molecule depends on its vibrational state \v). Therefore both, the bounding potential and 
'il>'^jj{R), depend on the vibrational quantum number v. The translational quantum numbers are 

= 0, 1, 2, . . . . The orbital angular momentum quantum number takes values L = A^, — 2, . . . , 1 
or 0, depending on the parity of A^, and the azimuthal quantum number is Ml = —L, —L+1, . . . ,L. 
The rotational wavefunctions, defined by the rotational quantum numbers J = 0, 1, . . . and Mj = 
— J, — J + l,...,+J, are given by the spherical harmonics YjMjid, 4>)- 

We use bipolar spherical harmonics with overall spherical rank A and component A/a 

Ff4(o,f),)= ci^itjMyLM.myjMA^s), (2) 

Ml,Mj 
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where C are the Clebsch-Gordan coefficients 6j]. Then the full wavefunction describing the motion 
of the H2 molecule is 

\vjNLKMt,) = r:\s)rNL{R)FkMS^, ^s) (3) 

where ^/^"''(s) = \v) is the vibrational wavefunction with the quantum number v. 

The Hamiltonian % for the trapped molecule includes coupling terms between vibrational, 
translational, and rotational motion. For simplicity, we neglect all matrix elements non-diagonal 
in V and introduce a parametric dependence on v, 

2 

= ''^"^-''°* + ^ + ^F(i?,0,l^s), (4) 
Im 

where ^^"^1^-''°* is the vibration-rotation Hamiltonian, p is the molecular momentum operator and 
m is the molecular mass of the diatomic molecule. = {v\ V{R, s, il, ilg) \v) is the potential energy 
of the hydrogen molecule in the vibrational state ipj/,^'^{s). The vibration-rotation Hamiltonian 
D^vib-rot -g (jiagonal in the basis set \vJNLAM\) with eigenvalues given by 

^^vib-rot ^ ^^^^ ^ ^/2) + BrJ{J + 1), (5) 

Br = Be — Oie{v + 1/2) — DeJ{J + 1), where ujq is the fundamental vibrational frequency, oce is the 
anharmonic correction and is the centrifugal correction to the rotational constant B^ 65, 6^ . 
We start from the general expansion of the potential in multipoles 

^V{R,n,ns)= Y.'^XmSR)FLS^^^s)- (6) 
and expand the radial part of potential ^V^^^ {R) in powers of R 

n 

where n > I and the parity of I and n is the same. 

y(R, s) is a scalar and transforms under fully symmetric representation Ag of the symmetry 
group Jft- The spherical harmonics A = 0,6,10,... transform like Ag of the symmetry group 
Ih [63]. We use fully spherical approximation of the potential, A = 0. Because A = and 
A = |Z — j|, |/ — j'l -|- 1, I + j it must be that / = j. 

The total potential is 

+ {XS'R' + XS'R') + {X^'R + X^'R') Fo^d + {XS'R' + %TR') Fo'l (8) 

if we limit our expansion to j = / = 2 and n = 4. The odd-j terms are not allowed by symmetry for 
H2 and D2 and thus the coefficients ^qq^ and "Vo(/'^ are zero for homonuclear diatomic molecules. 
If we set the constant off-set ^'o'^^oo ^ 

zero and write the perturbation part as 

= Xo''R'fSS + {Xo''R + Xo''R') + iXf'R' + %rR') Fo'l (9) 

and the isotropic harmonic term as ^V^ = '^^qq'^R'^Fqq, then the total Hamiltonian reads 

2m 
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TABLE I. Translation-rotation energies of a 3D spherical oscillator {vJNLAv\ U ~ ^T^^i^-rot \yjjs[LK) for a 
perturbation given by Eq.[9]and H by Eq.[lO]for few lower states. Here ''A"°'* = ^jj^^l^ , = 
''A^^* ^ and for ''/J see Eq.HH 

JAfLA (vJNLkv\ n - ^'^vib-rot |t,jAfLA) 



0000 /i^'i^g^ll + ^A"M] 
1001 + "A0''4] 

0111 7„^004] 

1110 ^"^([[f + + loC-A^aa + -A^a*)] 

nil f^^'^oll + t'"AOO* - 5 C'A^aa + -A^a*)] 

1112 ^l-l.^g^d + |"A004 + 1 {-A222 + ^.a224)] 

2111 ^"t^g^lf + |"A004 + 5 {"A222 + "A224)] 

2112 fi^^oii + 5 - 5 ("'^^^^ + ''^^^*)] 

2113 f^^^oll + 1"^°°* + f ("A222 + «A224)] 

0200 + 5"A004] 

0222 n"c.;J[| + f-A004] 

1201 ?i"tj^[- + S^A^"*] 

1221 + ^"A004 + 7C'A222 + |''A224)] 

1222 fl-^^oit + ^"A004 - 7C'A222 + |''A224)] 

1223 h^UJ^ll + I.-A004 + 2rA222 + |«A224)] 



The unperturbed Hamiltonian eigenvalues in the basis \vJNLAM\) are 

E'.JNLAM, = ''Ej^-'°' + ''^liN + 3/2), (11) 

where = ('^*Q'^/(27rm)) is the frequency for translational oscillations within the cavity. 

The meaning of different parts of the perturbation is explained by their influence on the energy 
levels of a harmonic 3D spherical oscillator, TablelH Translation-rotation coupling term ^qq'^ splits 
energy of \vJNLA) state into levels with different A, where A = \L — J\,\L — J\ + l, . . . ,L + J. For 
example, the N = L = J = 1 state is split into three levels with different total angular momentum 
A = 0,1,2. The ordering of levels depends on the sign of ^qq^- The anharmonic correction 
to translation-rotation coupling is ^qq^- If isotropic anharmonic correction 'T/qq'^ is positive the 
distance between energy levels increases with and this correction is different for the levels with 
same but different L. For example, for positive 'T/qc?^ the A^ = 2, L = level has higher energy 
than N = 2,L = 2 level. 

The length scale ^(3 = m'"ujQ /h (dimension m~^) of the radial part of a 3D spherical oscillator 
wavefunction is related to the expectation value of the center of mass amplitude in state |A^) as 

13 



(ATI \N) = ^r\N + 3/2) = hJ—^{N + 3/2). (12) 

V '^00 

Terms described by the translation-rotation coupling coefficients ^qq"^ and ^qq^ do not appear 
in TableU] because the first order correction to energies vanishes as the matrix element of Fqq is 
zero if diagonal in L or J. These terms mix states with different A^ and J. For example, the first 
excited rotational state J = 1, A^ = (expectation value of HD center of mass is on the cage 
center) has the state J = 0, A^ = 1 (expectation value of HD center of mass is off the cage center) 
mixed in [sj ]. The effect is that HD is forced to rotate about its geometric center instead of center 
of mass. 
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It was found by the 5D quantum mechanical calculation that the rotational quantum number 
J is almost a good quantum number for the homonuclear H2@C6o and D2@C6o and not for the 
heteronuclear HDQCeo ^4Q] . Indeed, ^qq^ mixes states with different J for homonuclear species 
as well but the effect is reduced compared to the effect of ^qo^'^- former case J it 2 and L it 2 

are mixed while in the latter case the J it 1 and L it 1 states that have a smaller energy separation 
are mixed. 

The states with different A are not mixed in the spherical approximation, i.e. the total angular 
momentum A = L + J is conserved and A is a good quantum number. The other consequence of 
the spherical symmetry is that the energy does not depend on Ma. Therefore it is practical to use 
a reduced basis and reduced matrix elements [g^ which are independent of Ma. This reduces the 
number of states by factor 2A + 1 for each A. 



Model Hamiltonian of H2@C7o 

A spherical approximation of the potential of a molecule trapped in C70 would be an oversimpli- 
fication because of the elongated shape of C70. The symmetry of C70 is D5/1, the distance between 
the centers of two capping pentagons {z direction) is 7.906 A . The diameter of the equatorial 
xy plane is 7.180 A[69|], similar to the diameter of the icosahedral sphere of Ceo, 7.113 A[70|]. The 
anisotropy of the potential of H2 inside C70 is supported by the 5D quantum mechanical calculation 



4ll ] what shows that the lowest translational excitation in the z direction is 54cm~^ and in the 
xy plane is 132 cm~^ while in Ceo it is 180 cm~^ and isotropic ^]. We derive from the IR spectra 
(see below) that the xy plane excitation energy is 151 cm~^, somewhat larger than theoretically 
predicted. 

Although the z axis translational energy in C70 is three times less than in the icosahedral Ceoj 
the effect of the C70 potential on the rotational motion is moderate. The splitting of the J = 1 
state is 7cm"^ what is relatively small compared to the rotational energy 120 cm~^ in this state 



41| 



To analyze the IR spectra of H2@C7o we use a simplified Hamiltonian where the translational 
energy is represented in the form of the sum of two oscillators, ID linear and 2D circular oscillator, 
and we do not consider anharmonic corrections and the translation-rotation coupling. 

The vibration-rotation energy 



.^vib-rot ^ ^ ^^(^ + 1 ) + J(J + 1) +^/t(3Jf - 2), (13) 



is the same as for H2@Ceo except the last term which accounts for the axial symmetry of the C70 
potential with the rotational anisotropy parameter [7l|]. For example, the three-fold degenerate 
J = 1 rotational state in 1^ symmetry is split in D^h symmetry and if > the = state is 
3^ below the twice degenerate Jz = ±1 rotational state. 

The translational part is added to the vibration-rotation Hamiltonian, Eq. ll3l and the total 
energy reads 



EljJ.^nln. = ''E^V"' + ^ '^^.(™ + 1) + ^ ''^^n. + -J- (14) 



Here the translational energy is written as a sum of two oscillators, a linear oscillator along the z 
axis with translational quantum h '"coj and a 2D (circular) oscillator 48|, |64 1 in the xy plane with 



translational quantum h^u^y. Quantum numbers and n are positive integers including zero and 
I = n,n — 2, . . . , —n + 2, — n. 

We will show below that the frequencies of z and xy translational modes, \i;J and ^uj'^y^ can 
be determined from the experimental data even though the translation-rotation coupling is not 
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TABLE II. Classification of energy levels of H2 inside the cage of C70 up to J = 1 and n = = 1 by 
irreducible representations Ti of the symmetry group D^h- 



J 


(nln^) 







(000) 




1 


(000) 


E[ 





(001) 




1 


(001) 







(110) 


K 


1 


(110) 


A\, A'^,E'^, E'{ 



known. We take the advantage of translation-rotation coupling being zero in the J = rotational 
state. The complication rises from the fact that the potential is different in the initial and final 
states of the IR transitions, f = and 1. However, this complication could be resolved if the energy 
of the fundamental vibrational transition f = — )• 1 (without change of n and n^) is known. 

The AJ = transition from the para-H2 ground state leads to two excitation energies in the 
IR spectrum, one for the z mode and second for the xy mode 

^?oo;OOi - ^ooo;000 = ^[^0 + + C^^y " + J (^o;^ - ^ivj)] , (15) 



^?oo;iio - ^o°oo;000 = + '^^y + C^j^y " + ^ ('^J " "o; J)] . (16) 



2 

1 

2 



Defining the fundamental para transition energy as -E'ioo;ooo ~ -^000000 ~ ^['^o + C^xy ~ ^^xy) + 
^ (^wj — "wj)] we may rewrite Eq.fTSl and [T6l as 

pO pO _ fci, ,T I pO pO /17^ 

-^100;001 -'^000;000 — "^2 -^100;000 -^000;000' K^' ) 

pO pO — f,\ ,^ _i_ po po 

-'^100;110 ~ -'^000;000 — " ^xy "r -^^lOOiOOO ~ -'^000;000- 

From these equations translational frequencies in the excited v = 1 state, and ^J^, can be 
determined without knowing the translation-rotation coupling. 

The classification of energy levels up to J = 1 and n = = 1 by irreducible representations 
Ti of the symmetry group Ds/j is given in Tablellll We get the irreducible representations Tj: 
L = ^ A'^ and L = 1 — t- A2 + E[ by subducting the translational states represented by spherical 
harmonics Y^Ml from the full rotational group 0(3) to the symmetry group D^h- ^'1 is the para- 
H2 ground state, n = I = Uz = 0. The first excited state of the z mode is Uz = 1 and A2. The 
first excited state of the xy mode n = Z = 1 is doubly degenerate E[. The full symmetry when 
translations and rotations are taken into account is Fj = Tj (8> T^'^\ For example, the ortho-il2 
ground state, J = 1 and (nlnz) = (000), is split into Jz = (Ag) and = ±1 {E[), see Tablellll 



Induced dipole moment of hydrogen in spherical environment 



IR light is not absorbed by vibrations and rotations of homonuclear diatomic molecules 65l |. 
IR activity of H2 is induced in the presence of intermolecular interactions, such as in the solid 
and liquid phases [tI, 73], in constrained environments [sl, [H^, 62, 74|, and in pressurized gases 
75l . [Taj. IR spectra of such systems are usually broad due to inhomogeneities in the system or 



due to random molecular collisions. As an exception, narrow lines are observed in semiconductor 
crystals 77 1 and solid hydrogen fl^ . 

An overview of collision-induced dipoles in gases and gas mixtures is given in the book by L. 
Frommhold [t^ . The confinement of the endohedral H2 introduces two differences as compared to 
H2 in the gas. First, the translational energy of H2 is quantized. In the gas phase it is a continuum 
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starting from zero energy. Second, the variation of the distance between H2 and the carbon atom 
is hmited to the translational amphtude of H2 in the confining potential. In the gas phase the 
distance varies from infinity to the minimal distance given by the collision radius. The selection 
rule AA^ = ±1 for the endohedral H2 follows from these two conditions, as shown below. 

Quantum mechanical calculations of induced dipoles are available for simple binary systems like 
H2-He, H2-Ar, and H2-H2. An extensive set of theoretical results for the H2-He system associated 
with the roto-translation electric dipole transitions, both in the vibrational ground state v = and 
accompanying the u = — )• 1 transition of the H2 molecule, can be found in 83-84|. Related to 



the fullerene studies are calculations of the dipole moment of COQCgo [63| and exohedral H2 in 
solid Ceo fsi]. 

We express the induced part of the dipole moment as an interaction between hydrogen molecule 
and Cgo- Another approach was used in |33l ] where the summation over 60 pair- wise induced dipole 
moments between H2 and carbon atoms was done. The relation between two sets of parameters 
was given 



We write the expansion of the dipole moment from the vibrational state v to v' in bipolar 
spherical harmonics and in power series of R as 

4,,(R,Q.) = ^E^'^'^t.^^^^.l^^.^^.)- (18) 

This is similar to the expansion of the potential discussed above, except the dipole moment is 
a polar vector while the potential is a scalar. The dipole moment transforms according to the 
irreducible representation Ti^ of the symmetry group Ih- The spherical harmonics of the order 



A = 1,5,7,... transform according to Tiu of the symmetry group Ih |67l |. We use A = 1 and 
are interested in ?; = — )• 1 transitions. In spherical symmetry it is sufficient to calculate one 
component of the dipole moment vector, m\ = 0, and if we drop the explicit dependence of d^,^ 
on V, v' and of ^ -^y^^ on A, ra\, the ra\ = component of the dipole moment reads 

do(R,J^.) = ^^^'^"i2"Fl^o(f^>f^^)- (19) 

As A = |Z — j|, | / — j| + 1, Z + j and A = 1 it must be that I = |j ± 1|. The possible combinations 
are {Ij) € {(01), (10), (12), (21), (23), . . .}. If we restrict the expansion up to n = / = 1, we get 



do(R, 0.) = + ^ {A'^'Fl^ + A'^'Fli) R. 



^010 

is zero for homonuclear diatomic molecules II2 and D2. It is the sum of the HD permanent 
rotational dipole moment in the gas phase and the induced rotational dipole moment inside Cgo- 
The selection rule is AiV = 0, AJ = 0, ±1, but J = is forbidden. 

The expansion coefficients A^^^ and A^'^^ are allowed by symmetry for both, homo- and hetero- 
nuclear diatomic molecules. The expansion coefficient A^^^ describes the induced dipole moment 
that is independent of the orientation of the diatomic molecular axis s, selection rule AA^ = 
0, ibl,AJ = 0, but = — )• is forbidden. A^'^^ describes the induced dipole moment that 
depends on the orientation of s, AA^ = 0, ±1, A J = 0, ±2, but A^ = — and J = — )• are both 
forbidden. All terms in Eq. [20] satisfy the selection rule AA = 0, ±1, but A = — t- is forbidden. 

Although HD has a permanent rotational dipole moment, the induced dipole moment dominates 
inside Ceo [34,]. 



IR absorption line intensity 331] is proportional to the thermal population of the initial state. 



If the thermal relaxation between para-B.2 and ortho-B.2 is very slow we can define a temperature 
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FIG. 2. Boltzmann population of the H2@C7o para ground translational state (000) (solid line) and first 
thermally excited (dotted line) state (001) and of the two ortho states, (dashed line) and E[ (dash-dot 
line), in the ground translational state (000), calculated using energy levels of the 5D quantum mechanical 
calculation '41]. Here the para states have J = and the ortho states have J = 1; irreducible representations 
are from Tablellll (Online version in Colour) 



independent fractional ortho and para populations of the total number of molecules M, where 
k = o,p selects ortho or para-Il2- Then the probability that the initial state \viJiNiLiAiM\.) is 
populated is 

where Ei is the energy of the initial state measured from the ground state v = N = and j runs 
over all para- (or ortho-) H2 states in the basis used, gj = 2Aj + 1 is the degeneracy of the energy 
level Ej. Please note that gi does not appear in the numerator because pi is the population of a 
individual state \viJiNiLikiMx^) although the reduced basis \vJNLA) is used. Hetero- nuclear HD 
has no para and ortho species and the coefficient in Eg. [20] must be set to one, = 1. 



Infrared absorption by H2@C7o 



The infrared absorption in H2@C7o is given by electric dipole operators do and d±i, which 
transform like A2 and E[ irreducible representations of D^h. At low T the symmetry-allowed 

transitions from the ground para state are A[ 4 A2 and A[ 4' E[. Two ortho states, A2 and E[, 

are populated at low T, Fig.[2l The transitions are A2 ^ A'^ and A2 —> E'l. The symmetry-allowed 

transitions from the E'^ state are E'^ E'^ and E'^ '^■^ A\, A'2, E2- 

The IR absorption line area is proportional to the Boltzmann population of the initial state, 
Eq.[20l The degeneracy gj of the energy level Ej is one or two in case of H2@C7o . Fig.[2]shows the 
Boltzmann population of the few first para and ortho levels in the v = vibrational state. The para 
state (000) population (solid line) starts to drop above 20 K as the first excited state (001) (black 
dashed line) at 54cm~^ above the ground state becomes populated and drains the population from 
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FIG. 3. Infrared absorption spectra of H2@C6o at 6 K are shown with the black hne. The green hne in panel 

(a) is the Gaussian fit of fundamental para and ortho transitions, AiV = A J ~ 0. The red line in p anels 

(b) , (c), (d) is the simulated spectrum with parameters taken from the model fit of 200 K spectra 33"]. The 
S{1)^ labeled green line in panel (c) is the Gaussian fit of the AA^ = 0, A J — 2 forbidden ortho transition, 
experimentally observed at 4612.5 cm^^ and predicted to be at 4613.1 cm^^ by our model [sl]. Line labels 
are the same as in Fig.|3^; the superscript to line label shows the final state A. 



the (000) state. The abrupt change of the population of ortho ground states above 2 K is because 
of the transfer of population from the non-degenerate A2 ground level (dotted line) to the doubly 
degenerate E[ (dash-dot line) 7cm"^ higher. The populations of A2 and E'^ states decrease above 
30 K as the first excited translational (001) ortho state get more populated. 



EXPERIMENT 



The endohedral complexes were prepared by "molecular surgery" as described in [111 . [10]; 
H2@C6o, D2@C6o and H2@C7o at Kyoto University, and HDOCgo at Kyoto University and 
Columbia University. The HDCiCgo sample was a mixture of the hydrogen isotopomers H2:HD:D2 
with the ratio 45:45:10. Since all Cqo cages are filled, the filling factor for HD is p = 0.45. The 
content of C70 sample was empty:H2:(H2)2=28:70:2 and the filling factor p = 0.7. Experimental 
absorption spectra were corrected for the filling factor. 

The para enriched sample was made at Columbia University using molecular oxygen as a spin 



catalyst for ortho-para conversion [26(]. Briefly, the H2@C6o adsorbed on the external surface of NaY 
zeolite was immersed in liquid oxygen at 77 K for 30 minutes, thereby converting the incarcerated 
H2 spin isomers to the equilibrium distribution at 77 K, ng/rip = 1. The liquid oxygen was pumped 
away and the endofullerene-NaY complex was brought back rapidly to room temperature. The 
para enriched H2@C6o was extracted from the zeolite with CS2 and the solvent was evaporated by 
argon. The powder sample in argon atmosphere and on dry ice arrived in Tallinn four days after 
the preparation. 

Powdered samples were pressed under vacuum into 3 mm diameter pellets for IR transmission 
measurements. Typical sample thickness was 0.3 mm. Two identical vacuum tight chambers with 
Mylar windows were employed in the IR measurements. The chambers were put inside an optical 
cold finger type cryostat with KBr windows. In the measurements, the chamber containing the 
pellet for analysis was filled with He exchange gas while the empty chamber served as a reference. 
Transmission spectra were obtained using a Bruker interferometer Vertex 80v with a halogen lamp 
and a HgCdTe or an InSb detector. The apodized resolution was typically 0.3 cm~^ or better. 

Transmission T^(u;) was measured as the light intensity transmitted by the sample divided by the 
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FIG. 4. Diagrams of the observed low temperature IR transitions in hydrogen isotopologs (a) H2, D2 and 
(b) HD incarcerated in Ceo- The initial states surrounded by a dashed box have the vibrational quantum 
number w = 0; all the final states are in the excited vibrational state t; = 1. Dashed lines in (a) are forbidden 
transitions {AN = 0) that are observed in H2@C6o but not in D2@C6o- The S(l) line (dotted) with AA^ = 1 
is not observed in D2@C6o- The line label Q{Ji) is used for AJ = 0, R{Ji) is used for AJ = +1 and S{Ji) 
for A J = +2 transitions where Ji is the initial state J. 



light intensity transmitted by the reference. The absorption coefficient a{u}) was calculated from 
the transmission Tj.{oj) through a{oj) = —d~^\n. [rr(a;)(l — with the reflection coefficient 

[(^ ~ l)/(^ + 1)]^ calculated assuming a frequency-independent index of refraction 85|, 77 = 2. 
Absorption spectra were cut into shorter pieces around groups of H2 lines and a baseline correction 
was performed before fitting the H2 lines with Gaussians. 



RESULTS AND DISCUSSION 



-^60 



The IR absorption spectra of H2@C6o at 6K together with the simulated spectra are shown in 
Fig. [3] and the diagram of energy levels involved in Fig. [3^. The simulated spectra are calculated 
using Hamiltonian, dipole moment and the ortho-para ratio parameters (see TablelTV]) obtained 
from the ffi of 200 K spectra [s^]- Temperature does not affect these parameters and the line 
intensities follow the Boltzmann population of initial states [s^]. Intensity of three lines, Q{Q)^ and 
Q{lY Fig-Ela) and S{1)^ in c) can not be simulated because the induced dipole moment theory 
does not describe AA^ = transitions. However, the position of these three lines was used to fit the 



Hamiltonian parameters |33l ]. The splitting of the ortho state N = J = 1 into A = 0, 1 and 2 by 
translation-rotation coupling is seen in Fig.[3]b). The para N = 1, J = 2 and ortho N = 1, J = 3 
states are split into three sublevels as well. However, because of the selection rule AA = ±1 only 
one para and one ortho S'-transition is IR active. 
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T = 6K 

njn^ = 2.890 ± 0.045 
T = 20K 

n/n = 1 .397 ± 0.074 
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FIG. 5. Spectra in the region of Q lines of a para enriched H2@C6o (dashed line) and non-enriched sample 
(solid line) at low T. The thermal equilibrium ortho-para ratio at 78 K where the para enrichment was done 
is rio/rip = 1. (Online version in Colour) 



The intensity of other lines in Fig.[3]b-d) is described accurately by two dipole moment parame- 
ters, A^'^^ and A^'^^, and ortho-para ratio Uo/np = 2.89 it 0.045, what is very close to the statistical 
value three, Tabl dlVl We confirmed the assignment of spectral lines to para- and ortho-Il2 by 
measuring the spectrum of a para enriched H2@C6o sample |3^. Low T spectra in the region of 
Q lines are shown in Fig.O The time delay between para enrichment and first IR measurement 
was 4 days. The 4255 cm^^ para line is stronger and other three ortho lines are weaker in the 
para-enriched sample, as compared to the non-enriched sample. 

The Fig. [5] deserves some attention. Lines Q(0)^ and Q{1)^ have a clear multicomponent struc- 
ture and it is different for the two samples. This is not possible in the approximation of spherical 
symmetry. Even not by considering the true icosahedral symmetry of Cgo because the lowest A 
value for the state which is split by the icosahedral symmetry is three while the initial and final 
states of the optical transitions under discussion have only A = 0, 1. Later measurements on the 
relaxed but initially para converted sample showed that the difference between the line shapes of 
the normal and para converted sample is not due to the different ortho-para ratio. Thus it is likely 
that the difference in the line splitting is due to a different impurity content of the two samples. 
However, it is not completely excluded that the crystal field or the distortion of the Cqq cage is 
responsible for the part of this splitting. Note that the Q{1)^ line of D2@C605 Fig.[6l bears similar 
splitting pattern as in H2@C6o- 

The D2@C6o spectrum, Fig.[6l is shifted to lower frequency compared to H2@C6o because of 
the heavier mass of D2. The spectrum has less lines than the H2@C6o spectrum. The missing 
transitions in D2@C6o are shown in red color in Fig.H^ and they belong to J-odd species which 
are minority for D2. 

The splitting of = J = 1 into A = 0, 1 and 2 states is similar in D2 and H2. The magnitude 
of the splitting, h'"uj'n ^ , see TablelH is less in D2 because \lIq is smaller although Vgo*^ and 

^Vqq'^ are similar for two isotopologs, Table lTVl Line (5(1)^ overlaps partially with (5(0)^ line, Fig. [HI 
However, this is not because the translation-rotation splitting is different for two isotopologs but 
because of a smaller anharmonic correction Oe of the rotational constant for D2 , Table lTVl 

The spectrum of HDOCgo is more simple, Fig. [71 because there are no para and ortho species 
and therefore only one state is populated at low T. There is one spectral line, not present in 
homonuclear dihydrogen, labeled R{0) in Fig. [Da and Fig. [71 HD has no inversion symmetry and 
the ban on A J = 1 transitions is lifted. The classification of this transition as A J = 1 is arbitrary 
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FIG. 6. Infrared absorption spectra of D2@C6o at 5K. The simulated spectrum is shown by red hue (para), 
blue line (ortho) and black thin line (sum of para and ortho spectra). The parameters for the simulated 
spectrum are taken from the model fit of the 90 K spectrum [34]. Line labels are the same as in Fig.Hh'; 
superscript to line label shows the final state A. 



first, because the weight of JNLA = 1201 is only 0.5 in this state [34i] and secondly, the change of 
translational state is AA^ = 2, not allowed by any of the dipole operator components in our model. 
The next component in weight 0.29 in the final state is JNLA = 2111 which makes this transition 
A^^^-active, A^ = 0— 7>1,J = 0— )-2 and therefore it would be more appropriate to classify it 
as an S line. Another interesting feature of the HD spectrum is the absence of AA^ = 0, A J = 1 
(0000 — )• 1001) transition although the induced dipole moment coefficient A^^^ gives a two orders of 
magnitude larger dipole moment than the permanent dipole moment of free HD, as was discussed 
in [3j]. This A J = 1 transition in HDQCgo is suppressed because there is an interference of two 
dipole terms, A^^^ and A^^^ which have opposite signs. The final state consists of 80% of the pure 
rotational state J = 1,A^ = and 18% of the the pure translational state = 1, J = 0. This 
mixed final state has matrix elements from the ground state for both A coefficients, A^^^ and A^^^, 
which nearly cancel each other. 

The observed low T IR transitions of H2, D2 and HDOCqo are collected in Table HIIl The 
content of the unperturbed \JNLA) state in the final state is above = 0.9 in homonuclear 
species in most cases while for HD it varies and could by as low as 0.5. The mixing of states 
is proportional to the energy separation of mixed states, Ej — Ei in the first order perturbation 
theory. Fqq couples states where Li = Lj and Jj = Jj. These states are far from each other as 
Li = Lj only if Nj = Ni ±2. The other term Fqq mixes Lj = Lj it 2 and Jj = Jj it 2 what are 
even further from each other for small N. It was found by Xu et al. [41] that J is almost a good 
quantum number in H2 and D2@C6o but not in HDOCqq. Indeed, the distance between states 
mixed decreases for HD as Fqq mixes Lj = Lj it 1 and Jj = Jj it 1 and for this Nj = Ni ±1. 

The vibrational frequency is redshifted from its gas phase value for H2 and D2, Table llVl 
The relative change of the frequency [wn{gas) — wo(C'6o)]/'^o(gas) depends on the cage and is 
independent of the hydrogen isotopomer [87|]. Based on our fit results, Table llVl we get that for 
the H2 a;o(C6o)/'^o(gas) = 0.9763 and for the HD a;o(C6o)/wo(gas) = 0.9773. We used the average 
of these two ratios to calculate the D2@C6o fundamental vibrational frequency, cjq = 2924 cm~^ 
what was necessary to fit the IR spectra to a model Hamiltonian [341] . At this point we cannot say 
how much of the redshift is caused by change in the zero-point vibrational energy and how much is 
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FIG. 7. Infrared absorption spectra of HDQCgo at 5K, black line. The parameters for the simulated 
spectrum (red line) are taken from the model fit of 90 K spectra 3J] . The broad lines indicated by arrows 
are not due to HD. Line labels are the same as in Fig.HlD; superscript to line label shows the final state A. 
(Online version in Colour) 



TABLE III. Experimental IR absorption line parameters of endohedral dihydrogen isotopologs in Cgo at 
5K. w is the transition frequency in wavenumber units, cm~^ and S is the integrated absorption line area 
in cm~^ units. JNLA are the quantum numbers of the main component with the weight \^,^\'^ in the 
final state v = 1. The of the main component JNLA in the initial v — state for H2 are 0.97 of 
0000 and 0.97 of 1001, for D2 are 0.99 of 0000 and 0.99 of 1001, and for HD is 0.95 of 0000. Integrated 
absorption cross-section a is used by some authors. The integrated line are S is transformed into a (unit 
cm/molecule) using a = SV/Af where Af is the number of hydrogen molecules in the volume V [HI. Here 
Af/V = 1.48 X 10^^ cm^"^, the number density of molecules in solid Cqq. * It was not possible to separate 
overlapping Q{Q)^ and Q{1)^ lines of D2 and the line area given is the sum of both transitions. 
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caused by the change of anharmonic corrections to the vibrational levels in the Cgo as our data set 
is limited to energy differences of v = and v = 1 levels only. Note that the vibrational frequency 
of H2@C6o is iOQ = 4062.4 cm~^, Table lTVl while the para vibrational transition is shifted up by 
9cm~^, to ujQ = 4071.4 cm"'^, Table lllll This is because the zero-point translational energies are 
different in the v = and v = 1 vibrational states [s^. 

The rotational constant Bg, the vibrational correction Oe, and the centrifugal correction Dg 
of the hydrogen inside Cgo and in free space are compared in Table lTVl The smaller than the gas 
phase value of Be may be interpreted as 0.8% (0.9%) stretching of the nucleus-nucleus distance s 
in H2@C6o (D2@C6o); as B(, ~ {l/s^). An attractive interaction between hydrogen atoms and the 
cage causes s to be longer. The elongation of the equilibrium proton-proton distance is consistent 
with the redshift of ujq |83] . However, the anharmonic vibrational correction to rotational constant, 
ae, is smaller inside the cage than in the gas phase. Here the cage has the opposite, repulsive, effect 



16 



TABLE IV. Values of the fitted parameters for HaQCeo 



60 



Tlie vibrational 



and for HD and DaQC, 

and rotational constants of a free molecule of the three hydrogen isotopologs are shown for comparison; gas 
phase Wo is calculated including all terms {v + 1/2)'^ up to fc = i 
for HDQCgo and D2@C6o- Parameter errors are given in (s^ . 



The parameter "Kpg^'* is set to zero 
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V = 
14.28 
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v = l 
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4.3 ■ 10 



10 



v = l 
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3.25 ■ 10-10 
7.1 ■ lO^" 



V = 
16.78 
2.16 ■ 10^1 
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HDOCeo 
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and reduces the elongation of the proton-proton distance in the excited v states when compared to 
H2 being in the free space. This is supported by the fact that within the error bars of D2@C6o 
is equal to of D2. The vibrational amplitude of D2 is less than of the H2 and therefore the 
repulsive effect of the cage becomes important at w > 1 for D2 . 

Among the rotational and vibrational constants of HDQCqo the centrifugal correction Df. to 
the rotational constant is anomalously different from its gas phase value compared to other two 
isotopomers, Table lTVl Positive Dg means that the faster the molecule rotates, the longer is the 
bond. We speculate that since the rotation center of HD inside the cage is further away from the 
deuteron, the centrifugal force on the deuteron increases and the bond is stretched more than in 
the free HD molecule. 

A similar system to the one studied here is exohedral H2 in Cgo- There H2 occupies the 
octahedral interstitial site in the Cqo crystal. The prominent features in the exohedral H2 IR 
spectra [s^ are the translational, AA^ = ±1, sidebands to the fundamental transitions, Av = 1 
and A J = 0, 2. The redshift of the fundamental vibrational frequency is about 60 cm~^, which is 
less than in H2@C6o where it is 98.8 cm~^. Also the separation of translational = and A^ = 1 
levels, approximately 120 cm~^, is smaller, when compared to 184.4 cm~^ in H2@C6o- It is likely 
that the main contribution to the latter difference comes from the larger van der Waals volume 
available for H2 in the octahedral site than in the Ceo cage. 



H2@C7o 

IR absorption spectra of H2@C7o measured at 5 K are shown in Fig. [8] and the scheme of energy 
levels with the low T transitions indicated with arrows is shown in Fig.[9j The transitions where 
A J = and the translational state changes by An = 1 or Ariz = 1 are labeled as Qxy{J) and Qz{J) 
where J is the rotational quantum number of the initial state. The transitions where A J = 2 are 
labeled by Sxy{J) and 5*2 (J). The z and xy translational modes have distinctly different energies. 
For the Q lines they are shown in different panels, Fig.[8]a and [8]b. The para and ortho S lines 
are well separated because of different rotational energies. The para S lines are shown in Fig.[8]c 
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FIG. 8. IR absorption spectra of H2@C7o at 5K. Features shown with arrows in panel (c) are hkely caused 
by impurities, not by endohedral H2. 




FIG. 9. Energy levels and observed IR transitions of H2@C7o from the ground para and ortho translational 
states in the w = state to the excited vibrational state v = \. The irreducible representations of the 
symmetry group D^h are given for few lower states and are listed in Tabljlll The ordering of the states is 



based on [4l|. ortho-Il2 transitions marked by a and b are Qz(l)'' and Qz(l)'' in TablelVl and start from 
the thermally excited state, which is E[ (000) according to [41|. There are three experimentally observed 
ortho-^2 transitions from the ground A2 (000) state: Qxy{^), Qxyi^Y, and Qxy{^Y- 



and ortho S lines are shown in Fig.[8]d. The Q lines cannot be sorted out into para and ortho 
that easily. An exception is the (5(1) transition, inset to Fig.[8]a. This is a pure vibrational 
transition, f = — )• 1, without the change of translational or rotational states. The corresponding 
(fundamental) para transition (5(0) would be at 4069.3 cm~^ according to the model, TablelVl We 
expect it to be much weaker than the ortho (5(1) transition, as in H2@C60) and for this reason it 
is not observed in the experiment. 
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FIG. 10. Normalized line area and Boltzmann population of para and ortho transitions of H2@C7o starting 
from the ground translational state, (a) para lines, (b) and (c) ortho lines. Solid line in (b) is the normalized 
intensity when the initial state is the ground state A2 {Jz ~ 0) and in (c) when the initial state is E'^ 
[Jz — ±1), 7cm~^ above Aj. This corresponds to 'fe > 0, Eq.[T3l Dashed line in (b) and (c) is the reversed 
situation when the ground state is doubly degenerate E'l and A2 is 7cm~^ above it, i.e. 'fc < 0. The 
normalized intensities are calculated using the translational energy levels from the 5D quantum mechanical 
calculation ^ij and Eq.[20l 



In Fig. [10] we plot the normalized line area and normalized Boltzmann population. Fig.llOla 
shows the T dependence of para lines. The assignment of the 4125.6 cm~^ line to Qz{^) and the 
4219.9 cm~^ line to Qxy{^) is supported by our model that is summarized in TablelVl and will be 
discussed below. 

The uniaxial symmetry of C70 splits the ortho ground state J = 1. This splitting is about 
7cm^^~ 10 K and creates a sharp T dependence of the ortho line intensity at low T as is shown in 
Fig.[2j This feature could be used to determine which of the Q hnes belong to the ortho species. 
Two lines, 4110.5 cm~^ and 4124.2 cm~^, have a T dependence consistent with the ortho transitions 
starting from the ground translational (000) state of E'^ symmetry, Fig.[TO]a. However, according 
to our model there is no transition close to 4124.2 cm~^ from the thermally excited ortho E[ state, 
TablelVl This might be due to our simplified model where the rotational anisotropy parameter 
is assumed equal in f = and v = 1 vibrational states. The rest of the ortho lines are shown in 
Fig.[TO]b. Several lines follow the same T dependence, the fundamental transition at 4063.2 cm~^, 
Qxyi'^) hnes 4213,4, 4224.7, and 4226.9 cm'^ and the 5^(1) line at 4678.8 cm'^. Lines Qz{l) at 
4116 and Sxy{l) at 4761 cm~^ have distinctly different T dependence than the rest of the lines. 
Such deviation could be due to the overlap of transitions starting from A2 and E[ symmetry states. 
However, none of the ortho lines behaves like the theoretical prediction, black line, so the possibility 
that the ground ortho state is doubly degenerate E[ and not the A2, i.e. °k < 0, is not ruled out 
by our data. The dashed black line shows the ground state population if the ground state is E[ 
and the A2 state is 7cm~^ above it. This T dependence is more close to the experimental situation 
than the solid black line, where > 0. The change of sign of affects intensities only a little if 
the transitions start from the thermally excited state, Fig.llOlc. where the red dashed line is the 
normalized Boltzmann population of A2 state i/this state is 7cm~^ above the ground state E[, 
i.e. if °K < 0. 

It is not unreasonable that < and the E[ (J^ = ±1) state is the ground rotational state 
of ortho-il2- It requires that there is an attraction instead of repulsion between H and C when H2 
is in the center of the C70 cage. The attraction is stronger for the H2 if its molecular axis is in 
the xy plane, Jz = ±1, where the distance between H an C is less than if the axis is along the z, 
Jz = 0, where the H-C distance is longer. Indeed, the attraction between C and H in H2@C6o was 
deduced from the redshift of vibrational energy and reduced rotational constant compared to H2 
in the gas phase ^]. However, it was found by 5D quantum mechanical calculation that A2 is the 
ground state 4l|]. More elaborate model with dipole moment parameters that describe the IR line 
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TABLE V. Low temperature IR transitions of H2@C7o observed at 5K from the ground vibrational state 
t; = to the excited vibrational state v = 1. J{nlnz) are the quantum numbers for the initial and final 
states and the labels indicate corresponding transitions in Fig. [8] and [9l The experimental spectra were 
fitted with Gaussians, Sexp is the experimental line area; the experimental line position Uexp could always 
be determined with a precision better than O.f cm~^. The model parameters, Ea lf 41 were determined 
as described in the text by matching the model frequency cj^od to the experimental frequency uj^xp for 
four lines indicated with error 0. The model parameters are vibrational frequency loq — 4069.1 cm~^, ID 
oscillator frequency along the z axis "ojj = 56.5 cm~^, 2D circular oscillator frequency '"ijo'^y — 150.9 cm~^, 
and the rotational anisotropy parameter 'fc = 3.1 cm^^. The rotational constant Bf. = 59.865 cm~^ and its 
corrections, — 2.974 cm~^ and = 0.04832 cm^-'^ were assumed to be the same as in H2@C6o- 



Frequency 

( cm^-'-) 


Sexp 

(cm-2) 




Initial 
V = 




Final 
v = l 


Label 


Error 
(cm-i) 




^mod 


J 


(nlriz) 


J 


{nlnz) 




4063.2 


4063.2 


0.146 ±0.005 


1 


(000) 


1 


(000) 


0(1) 







4069.1 







(000) 


1 


(000) 


Q(o) 




4110.3 


4110.3 


0.311 ± 0.014 


1 


(000) 


1 


(001) 


Qz(i)" 





4115.7 


4119.7 


0.824 ± 0.009 


1 


(000) 


1 


(001) 


Q.(i) 


-4.0 


4124.2 


4119.7 


0.045 ± 0.007 


1 


(000) 


1 


(001) 


Qz(i)'' 


4.5 


4125.6 


4125.6 


0.096 ± 0.007 





(000) 





(001) 


Q.(o) 





4213.3 


4214.1 


1.759 ±0.024 


1 


(000) 


1 


(110) 


Qxy(l) 


-0.8 


4220.0 


4220.0 


7.702 ±0.027 





(000) 





(110) 


Qxy(0) 





4224.2 


4214.1 


2.111 ± 0.066 


1 


(000) 


1 


(110) 


Qxy(lY 


10.1 


4227.2 


4214.1 


10.143 ± 0.067 


1 


(000) 


1 


(110) 


QxyW 


13.1 


4462.3 


4456.3 


0.570 ± 0.023 





(000) 


2 


(001) 


5,(0) 


6.0 


4543.4 


4550.7 


0.582 ±0.012 





(000) 


2 


(110) 


Sxy{0) 


-7.3 


4678.8 


4667.0 


0.831 ± 0.012 


1 


(000) 


3 


(001) 


5.(1) 


11.8 


4761.3 


4761.3 


1.900 ±0.019 


1 


(000) 


3 


(110) 


Sxy (1) 


0.0 



intensities together with the Boltzmann population could resolve the issue. 

TablelVl summarizes the assignment of transitions observed in the experiment. Since our model 
described in Theory sectioun did not include translation-rotation coupling we used the transitions 
where L and J are both zero or one of them is zero in the initial and final states. The fun- 
damental frequency ujq = 4069.1 cm^^ was chosen to match the experimental fundamental ortho 
transition (5(1); the rotational constant = 59.865 cm~^ and its corrections = 2.974 cm~^, 
Df. = 0.04832 cm^^ of H2@C6o were used. Next the two para transitions were matched, Qz{0) and 
Qxy{0) to obtain ^ioj = 56.5 cm~^, ^uj'^y = 150.9 cm^-"^. The 4220 cm~^ line was chosen as Qxy{0) 
because its T dependence, Fig.[lO]a), is similar to 6*2(0) and Sxy{0) line T dependence. The ortho 
ground state splitting is 9.4 cm~^ (7.4 cm~^ from the 5D calculation |4l{ |) with % = 3.1 cm~^ which 
was obtained by assuming that the (^^(l)" transition, (OOO)ii^^ — )• (001) A'^ in Fig.[9l is centered at 
4110.3 cm~^. It was assumed that the anharmonic corrections to translational energy and to % are 
same for the v = and v = 1 vibrational states. 

The variation of the experimental frequencies compared to the model frequencies is within 
iblOcm"^, TablelVl This is reasonable since our model did not include translation-rotation cou- 
pling and the translation-rotation coupling in H2@C6o creates splittings of this magnitude. The 
vibrational frequency ojq is larger in H2@C7o than in H2@C6o- This difference must be taken with 
a reservation because the correction due to the difference of zero point translational energies in two 
vibrational states is not taken into account for H2@C7o as it was done for H2@C6o [IS]- Again, a 
more elaborate model is needed for C70 than the two oscillator model described in Theory sectioUn 
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CONCLUSIONS 

IR absorption spectra of endohedral hydrogen isotopologs, H2 [sl, 33], D2 fsS], and HD 34 1 
in Ceo are informative, involving excitations of vibrations, rotations and translational motion of 
dihydrogen. The translational motion of the encapsulated molecule is quantized and coupled to 
its rotations because of the surrounding Cgo cage. The vibrational frequency of dihydrogen is 
redshifted compared to the gas phase value. Together with the smaller rotational constant it shows 
that the hydrogen bond is stretched inside the cage and there is an attraction between H (or D) 
and C. The heteronuclear HD does not rotate about its center of mass because of the surrounding 
cage. Different rotational and translational states are mixed and rotational quantum number J is 
not a good quantum number for HDQCeo- 

Our study shows that the vibrations and rotations of Cgo and the crystal field effects of solid 
Ceo are not important on the energy scale of IR measurements. If these effects are important their 
contribution to the IR spectra is the order of one wavenumber splitting of few absorption lines. 



Such small splitting is consistent with the NMR [281 . |29[ and heat capacity [301 ] results. 



C70 has ellipsoidal shape and this splits the translation-rotation states of H2. The transla- 
tional frequency is about 180 cm~^ in H2@C6o- In C70 this three dimensional mode is split into 
a two dimensional mode at 151 cm~^ and a one dimensional mode at 56cm~^. The 5D quantum 



mechanical calculation 411] is in very good agreement with this experimental result. 
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